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We study analytically the asymptotic late-time evolution of realistic rotating collapse. This is done by 
considering the asymptotic late-time solutions of Teukolsky’s master equation, which governs the evolution of 
gravitational, electromagnetic, neutrino and scalar perturbations fields on Kerr spacetimes. In accordance with 
the no-hair conjecture for rotating black holes we show that the asymptotic solutions develop inverse power- 
law tails at the asymptotic regions of timelike infinity, null infinity and along the black-hole outer horizon 
(where the power-law behavior is multiplied by an oscillatory term caused by the dragging of reference 
frames). The dumping exponents characterizing the asymptotic solutions at timelike infinity and along the 
black-hole outer horizon are independent of the spin parameter of the fields. However, the dumping exponents 
at future null infinity are spin dependent. The late-time tails at all the three asymptotic regions are spatially 
dependent on the spin parameter of the field. The rotational dragging of reference frames, caused by the 
rotation of the black hole (or star) leads to an active coupling of different multipoles. 
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I. INTRODUCTION 


The no-hair theorem, introduced by Wheeler in the early 
1970s [1], states that the external field of a black hole relaxes 
to a Kerr-Newman field characterized solely by the black- 
hole mass, charge and angular-momentum. The mechanism 
for this relaxation process of neutral fields was first analyzed 
by Price [2] for a nearly spherical collapse. The physical 
mechanism by which a charged black hole, which is formed 
during a gravitational collapse of a charged matter, dynami- 
cally sheds its charged hair was first studied in [3,4]. How- 
ever, these analyses were restricted to spherically symmetric 
backgrounds, i.e. to the Schwarschild and _ Reissner- 
Nordstrom black holes. On the other hand, the physical pro- 
cess of stellar core collapse to form a black hole is expected 
to be non-spherical in nature because of stellar rotation. The 
analogous problem of the dynamics of massless perturba- 
tions outside realistic, rotating black holes is much more 
complicated due to the lack of spherical symmetry. 

Recently, this problem was addressed numerically [5,6] 
and analytically for the case of scalar perturbations [7]. How- 
ever, the evolution of higher-spin fields (and in particular 
gravitational perturbations) outside realistic, rotating black 
holes was not analyzed analytically so far. In addition, the 
works done so far were restricted to the asymptotic regions 
of timelike infinity [5-7] and along the black-hole outer ho- 
rizon [7]. The asymptotics at null infinity of rotating collapse 
was not studied so far. In this paper we give for the first time 
a full analytic analysis of the late-time evolution of realistic 
rotating collapse. Our analysis considers massless fields with 
arbitrary spin (and in particular gravitational perturbations). 
In addition, we consider the late-time evolution at all the 
three asymptotic regions: timelike infinity, future null infin- 
ity and along the black-hole outer horizon. 

There are two different approaches to the study of pertur- 
bations of Kerr spacetimes. The first is to consider perturba- 
tions of the metric functions. However, this direct approach 
leads to gauge-dependent formulations. An alternative ap- 
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proach is to consider curvature perturbations (perturbations 
of the Weyl scalars). Based on the tetrad formalism by New- 
man and Penrose, Teukolsky derived a master equation gov- 
erning the perturbations of Kerr spacetimes [8,9]. The main 
goal of this paper is to analyze the asymptotic late-time evo- 
lution of a realistic rotating collapse. This is done by study- 
ing the asymptotic late-time solutions of the Teukolsky equa- 
tion. 

The plan of the paper is as follows. In Sec. II we give a 
short description of the physical system and formulate the 
evolution equation considered. In Sec. III we formulate the 
problem in terms of the black-hole Green’s function using 
the technique of spectral decomposition. In Sec. IV we study 
the late-time evolution of perturbations of Kerr spacetimes, 
i.e., the asymptotic late-time solutions of the Teukolsky 
equation. We conclude in Sec. V with a brief summary of 
our results. 


II. DESCRIPTION OF THE SYSTEM 


We consider the evolution of gravitational, electromag- 
netic, neutrino and scalar (massless) perturbations fields out- 
side a rotating collapsing star. The external gravitational 
field of a rotating object of mass M and angular-momentum 
per unit-mass a is given by the Kerr metric, which in Boyer- 
Lindquist coordinates takes the form 


5 2Mr\_, 4Marsin? 0 
ds“=| 1 5 t+ 5 tdo 
5 2Ma’rsin? 0 
-p4 Sd -sin 9 rta’ +— o dg’, 
(1) 


where > =r?+a7cos?@ and A=r?—2Mr+a?. Throughout 
this paper we use G=c=1. 
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Using the Newman-Penrose formalism, Teukolsky [8,9] 
derived a master equation that governs the evolution of per- 
turbations of the Kerr spacetime 


r’ +a’)? 6? 4Mar 6 
See —a’sin’ ð ve { 
A ðt A dtdg 
2 2 
a 1 joe eee meee! 
A sin? @|dg7 or or 
1 ð|. gee 
sind 00\""". 40 
a(r—m)  icos0|dw M(r?—a’) 
2s a) S r 
A sin | ðe A 
op 
—iacosé ae 
+ (s*cot?6—s) p=0, (2) 


where the parameter s is the spin-weight of the field. The 
Teukolsky’s master equation is valid for scalar perturbations 
fields (s=0), neutrino perturbations fields (s = + 1/2), elec- 
tromagnetic perturbations fields (s=+1) and gravitational 
perturbations (s= +2). The field quantities y% which satisfy 
this equation (for the various values of the spin parameter s) 
are given in [9]. Resolving the field in the form 


o0 


pHa +a)" > | Prei”? (3) 


m=—% 


one obtains a wave-equation for each value of m: 
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where the tortoise radial coordinate y is defined by dy 
=[(r*+a’)/A]dr. The coefficients B;(r, 0) are given by 


Aa’sin? 0 
Bise ae ge (5) 
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and 
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B3(r,0)=;\2(s+ 1)(r—M)[s(r—M)A7! 


a? 1 
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+(steoto-9)| (7) 


(r?+a7)?" 
Ill. FORMALISM 


The time evolution of a wave field described by Eq. (4) is 
given by 


We) =20[ [Berteeesn Veo 


Gare see te 40) 
+B,(z')G(z,z'st)¥(z',0)}sind'd6'dy’, (8) 
for t>0, where z stands for (y,0). The (retarded) Green’s 
function G(z,z’;t) is defined by 
2 2 
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= d(t) y—y") (9) 


The causality condition gives us the initial condition 
G(z,z';t)=0 for t<O. In order to find G(z,z’;t) we use the 
Fourier transform 


Gi(y.y's0),SH(0.aw)=27 | [Tezis Oaw) 
0 JO 


X sinde!'d dt, (10) 


where „S7'(0,aw) are the spin-weighted spheroidal harmon- 
ics which are solutions to the angular equation [9] 


1 əf. ð pags dg an ine 
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For the aw=0 case, the eigenfunctions ,5S/'(@,aw) reduce 
to the spin-weighted spherical harmonics ‚Y7 (0, ¢%) 
= S7 (0)e™?, and the separation constants ,Aj"(aw) are 
simply ,A”=(1-s)(1+s+1) [10]. 

The Fourier transform is analytic in the upper half w- 
plane and it satisfies the equation [9] 
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d? |K*—2is(r—M)K+A(4irws—h) 
dy?" (Cray 


2 dH rand f 
-H E G(y,y";w) 


=ô(y-y'), (12) 


where K=(r?°+a°)w—am, N=A+a°w?—2amw and H 
=s(r—M)/(r?+a°)+rA/(r?+a®)?. G(z,z';t) itself is 
given by the inversion formula 


1 œ +ic X 
Glz) = A l G(y.y';w) 
=f A 


tj 


x S7 (0,aw) S7 (0',aw)e™ dw, (13) 
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s|). 


Next, we define two auxiliary functions Wi(z.w) and 


where c is some positive constant and lọ =max(|m 


’ 


F,(z,w) which are (linearly independent) solutions to the 
homogeneous equation 


d? ii K?—2is(r—M)K+A(4irws— X) 
dy* Geta 
dH 


eye aS 
dy 


[towo i=1,2. (14) 


The two basic solutions that are required in order to build the 
black-hole Green’s function are defined by their asymptotic 
behavior: 


ATS eT iY y>- 
E E T y>, (15) 
and 
Bouw) Ae Bp (wA, y>, 
P(y,w)~ y M= siwy yoo (16) 


where k=w-mw,,w,=a/(2Mr,). Let the Wronskian be 


Ww)=W(¥,,¥,)=P,¥,,-P.¥ (17) 


Ly? 


where W(w) is y-independent. Thus, using the two solutions 


W, and W,, the black-hole Green’s function can be ex- 
pressed as 


1 VP i(y,w)Po(y',w), ysy’, 


Ww) Wily’ w)Vo(y,w), yoy’. 
18) 


In order to calculate G(z,z’;t) using Eq. (13), one may 
close the contour of integration into the lower half of the 
complex frequency plane. Then, one finds three distinct con- 
tributions to G(z,z’;t) [11]: 

(1) Prompt contribution. This comes from the integral 
along the large semi-circle. This term contributes to the 
short-time response. 

(2) Quasinormal modes. This comes from the distinct sin- 


Giy.y'sw)= 


gularities of G(y,y';w) in the lower half of the complex 
w-plane, and it is just the sum of the residues at the poles of 
G(y,y';w). Since each mode has Imw <0 this term decays 
exponentially with time. 

(3) Tail contribution. The late-time tail is associated with 
the existence of a branch cut (in V,) [11], usually placed 
along the negative imaginary w-axis. This tail arises from the 
integral of G(y,y';w) around the branch cut and is denoted 


by GÏ. As will be shown, the contribution G© leads to an 
inverse power-law behavior (multiplied by a periodic term 
along the black-hole outer horizon) of the field. Thus, GF 
dominates the late-time behavior of the field. 

The present paper focuses on the late-time asymptotic so- 
lutions of the Teukolsky equation. Thus, the purpose of this 
paper is to evaluate G°(z,z';1). 


IV. LATE-TIME BEHAVIOR OF REALISTIC 
ROTATING COLLAPSE 


A. The large-r (low-w) approximation 


It is well known that the late-time behavior of massless 
perturbations fields is determined by the backscattering from 
asymptotically far regions [12,2]. Thus, the late-time behav- 
ior is dominated by the /ow-frequencies contribution to the 
Green’s function, for only low frequencies will be backscat- 
tered by the small effective potential (for r>M) in Eq. (14). 
Thus, as long as the observer is situated far away from the 
black hole and the initial data has a considerable support 
only far away from the black hole, a large-r (or equivalently, 
a low-w) approximation is sufficient in order to study the 
asymptotic late-time behavior of the field [13]. Expanding 
Eq. (14) for large r one obtains [using \=(/—s)(i+s+1) 
+ O(aw) and neglecting terms of order O(w/ r°) and higher] 


(1+1) 


d? 4Mw?+2isw 
+ -z P=0. (19) 
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It should be noted that this equation is the correct form of 
Eq. (66) of Ref. [13] (the term s?/r? should not appear in that 
equation). 

We introduce a second auxiliary field ¢ defined by 


Paper d(x); (20) 


where x=—2iwr.  (z) satisfies the confluent hypergeo- 


metric equation 


d? d 
ratta) (141-2809) |80. 
(21) 


where 


is 
a=M+—. (22) 
2w 
Thus, the two basic solutions required in order to build the 
Green’s function are (for r>M, al) 


W,=Ar't eM (14+5+1-—2iwM,21+2,-2iwr), 
(23) 


and 


W,=Br'* e™'U(1+s5+1—2iwM,214+2,-2iwr), 
(24) 


where A and B are normalization constants. M(a,b,z) and 
U(a,b,z) are the two standard solutions to the confluent hy- 
pergeometric equation [14]. From these solutions it is clear 
that the black-hole rotation parameter a is irrelevant in the 
context of the asymptotic form of the Green’s function (to 
leading order in the inverse time). However, as will be 
shown below, the rotational dragging of reference frames, 
caused by the rotation of the black hole (or star), has (two) 
important effects (which are absent in the non-rotating case) 
on the asymptotic late-time evolution of rotating collapse. 
Hence, in order to find the asymptotic form of the Green’s 
function it is sufficient to analyze the late-time solutions in 
the spherically symmetric limit (i.e., on the Schwarschild 
background). 

A note is needed here on the relations between the various 
equations governing the perturbations of the Schwarschild 
spacetime. There are two different approaches to the study of 
linearized gravitational perturbations: by considering metric 
perturbations, or, alternatively, by considering perturbations 
of the Weyl scalars. Equations governing metric perturba- 
tions of the Schwarschild spacetime were derived by Regge 
and Wheeler [15] (for odd-parity perturbations) and by Zer- 
illi [16] (for polar perturbations). An alternative approach, 
based on the tetrad formalism by Newman and Penrose was 
used by Teukolsky [7,8] to derive a separable wave equation 
for perturbations of the Weyl Scalars. These scalars, con- 
structed from the Weyl tensor in a given tetrad basis, char- 
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acterize the gravitational field in vacuum. They allow a more 
convenient approach to the study of gravitational perturba- 
tions, due to their scalar nature. As was first shown by Chan- 
drasekhar [17], the a— 0 limit of the Teukolsky equation, the 
Bardeen-Press equation [18], can easily be transformed into 
the Regge-Wheeler equation. However, this equivalence is 
restricted to the frequency domain. In other words, any solu- 
tion (in the frequency domain) of the Bardeen-Press equation 
can be transformed in a trivial manner into a solution of the 
Regge-Wheeler equation, and vice versa. However, these so- 
lutions have, of coarse, a different w-dependence [as can be 
verified for example from their asymptotic forms Eqs. (23) 
and (24)]. Hence, this difference in the w-dependence of the 
solutions, when integrated in the complex frequency plane 
(in order to obtain the temporal dependence of the solutions) 
may lead to a non-trivial difference in the time dependence 
of the solutions. Thus, it is not a priori guaranteed that the 
asymptotic late-time solutions of Teukolsky’s equation and 
the Regge-Wheeler equation (studied by price [2] for the 
spherically symmetric case) will have the same temporal be- 
havior. (Indeed, it is found that the dumping exponents, de- 
scribing the fall-off of the asymptotic solutions at future null 
infinity are spin-dependent, a phenomena not found for the 
Regge-Wheeler asymptotic solutions; see Sec. IV C below.) 
Furthermore, one can easily verify that even the a—0 limit 
of the Teukolsky equation depends on the spin parameter s 
[see Eq. (2) and its asymptotic form Eq. (19)]. Hence, it is 
not surprising that the asymptotic late-time solutions of the 
Teukolsky equation will turn out to be spin-dependent. 

The function U(a,b,z) is a many-valued function, i.e., 
there will be a cut in V,. Using Eq. (13), one finds that the 
branch cut contribution to the Green’s function is given by 


1 f-i 
C R E ' 
Gead ae, J, io" 


Wiy,wer™) Wa(y,w) 
W(we27") W(w) 


x S(0,aw) Si(0',aw)e dw. (25) 


(For simplicity we assume that the initial data has a consid- 
erable support only for r-values which are smaller than the 
observer’s location. This, of course, does not change the 
asymptotic late-time behavior.) 

Using the fact that M(a,b,z) is a single-valued function 
and Eq. 13.1.10 of [14] (taking the b—k limit, where k is an 
integer), one finds 


P (rwee™)=V (rw), (26) 


and 
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"o ES B —1)?!2 ri 
Pa(r,we?™)=WVa(r,w) + = 


a OG DI CIs- 2iwM) OY) (27) 


Using Eqs. (26) and (27) it is easy to see that 
W(we?™)=W(w). (28) 
Thus, using Eqs. (26), (27) and (28), we obtain the relation 


W5(r,we?™) Vi(rw) B (—1)7/2 wi WV (r,w) 
Wwe") Ww) sO «(274+ 1)!T (14+ s—2iwM) Ww) ` 2) 


Since W(w) is r-independent, we may use the large-r asymptotic expansions of the confluent hypergeometric functions (given 
by Eqs. 13.5.1 and 13.5.2 in [14]) in order to evaluate it. One finds 


AB(-1)!'*1(2141)!w7 771 


WW)=i F541 —2iwM) 22 (30) 


(Of course, using the |z|—0 limit of the confluent hypergeometric functions, we obtain the same result.) Thus, substituting 
Eqs. (29) and (30) in Eq. (25) one finds 


= o (=1)+!2% mo o l 
C hig 1 1 21+1,-iwt 
Gln 2 Samael) TOI 2iwM) YOO EO wS 8,4) Si aw wr! e idw. 


(31) 


As was explained, the late-time behavior of the field should follow from the low-frequency contribution to the Green’s 
function. It is clear that all the rotation-dependent terms of Eq. (11) can be neglected in the aw—0 limit. Hence, one may 
replace the functions ,S/"(6,aw) by the spin-weighted spherical harmonics ,Y/"(6,@). In addition, we use the approximation 


Td+s+1-—2iwM) 
T'(-Il+s—-2iwM) 


=2iwM(—-1)!°*!(1+s)!(I—-s)!, (32) 


which is valid for w—0. With these substitutions Eq. (31) becomes 


og ex EM(- 17-5274 1048) 11-8)! oo ie SoN 
G(z,2 = 2 TATOI+ DIE sY (0) V7 (0 f Vi (yw) Pi (y ,w)w” te "dw. (33) 


B. Asymptotic behavior at timelike infinity 


First, we consider the asymptotic behavior of the fields at timelike infinity i, . As was explained, the late-time behavior of 
the field should follow from the low-frequency contribution to the Green’s function. Actually, it is easy to verify that the 
effective contribution to the integral in Eq. (33) should come from |w|=O(1/t). Thus, in order to obtain the asymptotic 


behavior of the field at timelike infinity (where y,y'<t), we may use the |w|r<1 limit of V,(r,w). Using Eq. 13.5.5 from 
[14] one finds 


Wi (rw)=Ar'tl, (34) 


Substituting this in Eq. (33) we obtain 


=e Z iM(—1)?!-82744 (1+ 5)!—s)! 
Caz! = 2 "AESI 


SYO YE Cy j w”! +2e~iWidw, (35) 
0 


Performing the integration in Eq. (35), one finds 


o 


; M(-1)!'-8*!27/*1(7+45)'(1—s)1(21+2)! 
ezi ait DIP 
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C. Asymptotic behavior at future null infinity 


Next, we go on to consider the behavior of the fields at future null infinity scri, . It is easy to verify that for this case the 
effective frequencies contributing to the integral in Eq. (33) are of order O(1/u). Thus, for y—y'’<t<2y—y' one may use 


the |w|y’<1 asymptotic limit for V,(y’,w) and the |w|y>1 (Imw<0) asymptotic limit of V,(y,w). Thus, 
Fiy, w)=Ay "t, (37) 


and 


eT SES AWM Jyp) Se PEM gee rer 


Vi (y,w)=Ae (2141)! (38) 


T(d-s+1+2iwM) > 
where we have used Eqs. 13.5.5 and 13.5.1 of [14], respectively. Integrating Eq. (33) with the aid of Eqs. (37) and (38) one 
finds (for v>u) 


2 M(-1) 5+ 124 4s)1—84+2)1(2142) g Lag 
OO ry,? —s,—(I-s+2) 
£ m(1+2)! sY O) YFCO')y" v “u . (39) 


D. Asymptotic behavior along the black-hole outer horizon 


Finally, we consider the behavior of the fields at the black-hole outer-horizon r, . While Eqs. (23) and (24) are (approxi- 
mated) solutions to the wave equation (14) in the r>M,|a| case, they do not represent the solution near the horizon. As y 
—— œ the wave equation (14) can be approximated by the equation [9] 


2is(r,-M)k  s*(r,—M)? 
2 = 
Foyt s 2Mr, (2Mr.)? V=0. (40) 
Thus, we take 
E (y,w)=C(w) Ae, (41) 


and we use Eq. (37) for V,(y’,w). In order to match the y<—M solution with the y>M solution we assume that the two 
solutions have the same temporal dependence (this assumption has been proven to be very successful for neutral [19] and 
charged [3] perturbations on spherically symmetric backgrounds). In other words we take C(w) to be w-independent. In this 


case one should replace the roles of y’ and y in Eq. (33). Using Eq. (33), we obtain 


o0 


M(~-1)!~8*!27"*1(74+s)1(1—s)1(214+2)! 


<Y 0) YF ( O')A Sy ttt Ipimwiyy =(21+3), (42) 
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where Tọ is a constant. 


E. Coupling of different multipoles 


The time evolution of a wave field described by Eq. (4) is 
given by Eq. (8). The coefficients B,(r,0) and B3(r,0) ap- 
pearing in Eq. (4) depend explicitly on the angular variable 0 
through the rotation of the black hole (no such dependence 
exist in the a=0 case). This angular dependence of the co- 
efficients leads to an active interaction between different 
multipoles (characterized by different values of /). This cou- 
pling between different multipoles is physically caused by 
the rotational dragging of reference frames (due to the rota- 
tion of the black hole). Thus, even if the initial data is char- 
acterized by a pure multipole / (i.e. by a certain spin- 
weighted spherical harmonic ,Y,) other multipoles would be 
generated dynamically during the evolution. In other words, 
the late-time evolution of a spin-s field is dominated by the 


lowest allowed multipole, i.e., by the /=|s| multipole, re- 
gardless of the angular dependence of the initial-data (pro- 
vided that the initial-data contains modes with |m|<|s|). 
This phenomena of coupling between different multipoles 
has been observed in numerical solutions of the Teukolsky 
equation [5,6] (a similar phenomena is known in the case of 
rotating stars [20]). 


V. SUMMARY 


We have studied the asymptotic late-time evolution of 
realistic rotating collapse. This was done by considering the 
asymptotic late-time solutions of Teukolsky’s master equa- 
tion which governs the evolution of gravitational, electro- 
magnetic, neutrino and scalar (massless) perturbations fields 
on Kerr spacetimes. 

Following the no-hair conjecture for rotating black-holes 
we have shown that the asymptotic solutions develop inverse 
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power-law tails at timelike infinity, at null infinity and along 
the black-hole outer horizon (where the power-law behavior 
is multiplied by an oscillatory term, caused by the dragging 
of reference frames at the event horizon). 

The dumping exponents, describing the fall-off of the 
asymptotic solutions of Teukolsky’s equation at timelike in- 
finity and along the black-hole outer horizon are independent 
of the spin parameter of the field. However, we have shown 
that the dumping exponents at future null infinity are spin- 
dependent. Moreover, the asymptotic late-time solutions of 
Teukolsky’s equation at all the three asymptotic regions have 
a spatial dependence on the spin parameter of the field. 

The dumping exponents are independent of the rotation 
parameter a of the black hole (or star). However, the rota- 
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tional dragging of reference frames, caused by the rotation of 
the black hole has two important effects (which are not 
found in the nonrotating case) on the asymptotic late-time 
evolution of rotating collapse: The power-law tail along the 
black-hole outer horizon is multiplied by an oscillatory term. 
In addition, the rotation leads to an active coupling of differ- 
ent multipoles. Hence, the late-time evolution of a spin-s 
field is dominated by the lowest allowed multipole, i.e., by 
the /=|s| multipole. 
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